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We present the first results of a high-resolution study of the distribution and kinematics of molecular gas in two nearby 
radio galaxies, 4C 31.04 and 3C 293, representative of two different stages of evolution in radio-loud active galactic nuclei 
(AGN). These observations, conducted with the IRAM Plateau de Bure Interferometer (PdBI), map with unprecedented 
spatial resolution (^0.5^^-1^0 and sensitivity the emission and absorption of key molecular species such as CO, HCN and 
HCO^. We report on the detection of a kinematically disturbed and massive (M^as ^lO^^M©) molecular/dusty disk of 
^1.4 kpc-size fueling the central engine of the compact symmetric object (CSO) 4C 31.04. We also report on the detection 
of a massive (Mgas ^lO^^M©) regularly rotating ^7 kpc-size disk in the PR II radio galaxy 3C 293. A complex system 
of molecular line absorptions is detected against the mm-continuum source of this galaxy (AGN and jet). We compare the 
properties of the molecular disks in the two sources and discuss them in the light of the different theories describing the 
evolution of radio galaxies. 

© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



1 Introduction 

Molecular gas dominates the total mass budget of neutral 
gas in the circumnuclear regions of spiral galaxies. Molecu- 
lar line observations are thus required to provide a complete 
view of the distribution and kinematics of neutral gas, a key 
to understanding how and for how long activity can be sus- 
tained in active galactic nuclei (AGN). Furthermore, AGN 
are able to inject vast amounts of energy in the surround- 
ing Interstellar Medium (ISM) of their hosts, most partic- 
ularly, in their circumnuclear molecular disks. Therefore, a 
complete understanding of the complex status of ISM re- 
quires also mm-line observations to probe the feedback in- 
fluence of activity on the energy balance/redistribution of 
the ISM in AGN. The central gas disks of radio-loud AGN 
(hosted by early-type galaxies) may well be also predomi- 
nantly molecular. Unfortunately, most of the molecular line 
observations of radio galaxies to date have been plagued by 
insufficient sensitivity or low spatial resolution (e.g., Evans 
et al. 2005; Prandoni et al. 2007; SaripalH & Mack 2007). 
As molecular lines are expected to be in emission but also 
in absorption (against the strong radio continuum sources) 
in radio galaxies, high spatial resolution and high sensitivity 
observations are needed to, at the least, reduce the bias due 
to the confusion of absorption and emission in the beam. 
Moreover, the combination of emission and absorption lines 
is paramount to get a complete 2D picture of the controver- 
sial dynamical status of the gas disks in radio galaxies. Con- 
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elusions based only on absorption line studies can be very 
misleading if we want to establish if the gas kinematics are 
dominated by inflow motions, outflow motions or circular 
rotation. 

In this paper we present the first results of a molecular 
line study of two nearby radio galaxies: the compact sym- 
metric object (CSO) 4C 3 1 .04 and the Fanaroff-Riley type II 
object (FR II) 3C 293. The study is based on high-resolution 
(0.5''-1'0 observations done with the Plateau de Bure inter- 
ferometer (PdBI) that make use of several species (like CO, 
HCO+ and HCN) adapted to probe the emission and the 
absorption of molecular gas in these sources. A study of the 
content, distribution, and kinematics of cold gas via molecu- 
lar lines is a key to elucidating the evolutionary link between 
young radio galaxies, thought to be represented by CSO, 
and the more evolved radio galaxy class, of which FR I and 
FR II are thought to be fair representatives. 



2 Young radio galaxies: the case of 4C 31.04 

4C 31.04 is a nearby CSO (z--0.06) lying at the nucleus of 
the giant elliptical galaxy MCG 5-8-18. HI absorption was 
first detected in this CSO by Van Gorkom et al. (1989) (see 
also Conway 1996). The HST maps of Perlman et al. (2001) 
reveal that the optical image of 4C 31.04 is permeated with 
dusty features; these consist of a disk and a spiral structure. 
The central disk is perpendicular to the axis of the ~100 pc- 
size radio source identified in the 5 GHz VLBA image of 
Giovannini et al. (2001). 
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Fig. 1 a) The Imm-continuum PdBI map of 4C 31.04 (contours: 3cr, 4a, 7a to 55a in steps of 6cr; l<j=0.8 mJy beam~^) 
is overlaid on the R-H color image of Perlman et al. (2001) from HST (color scale). (Aa, A(^)-offsets are relative to the 
AGN locus, b) The Imm-continuum map, degraded to the resolution of the 3mm observations (contours: 10% to 90% 
in steps of 10% of the peak value), is overlaid on the spectral index map (in color scale) of the continuum emission 
derived between 84 GHz and 218 GHz. c) The HCO+ (1-0) line maps of 4C 31.04 for the three velocity channels de- 
fined in Fig.[2]3: I(emission)(black thin contours; with maximum at B), II(absorption)(white contours; with maximum 
at the AGN), III(emission)(black thick contours; with maximum at A). Levels are 3.5a, 4a, and 4.5a for channel I 
(lcr(I)=0.06Jykms"^), -40cr, to -lOOa in steps of -lOa for channel II (lcr(II)=0.05Jy kms"^) and 5a to 13cr in steps 
of 2a for channel III (lcr(III)=0.08Jy kms~^). HCO+ contours are overlaid on the F702W HST image of Perlman et 
al. (2001) (color scale in mag.pixel"^). Ellipses represent the beams at 218 GHz(a)) and 84 GHz(b)). See GB07 for details. 
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Fig. 2 a-b) Spectrum of the HCO+(1-0) line emission and absorption at the central offset (AGN) of 4C 31.04. A strong 
absorption line is detected against the 3mm continuum source. Velocity scale is relative to z=0.0592. b) A zoom on the 
HCO+(1-0) spectrum of a). Intervals I-to-III and Vsys (2^=0.0602) value are highlighted, c) A high-resolution version of 
the HCO+(1-0) absorption profile of a); the three distinct velocity components identified in the spectrum are all blueshifted 
with respect to Vsys- See GB07 for details. 



Garcia-Burillo et al. (2007) (hereafter GB07) have used 
the PdBI to probe the gas and dust content of 4C 31.04 by 
observing the HCO+(1-0) and ^^00(2-1) molecular lines 
and their underlying continuum emission with high spatial 
resolution (<1'0- The continuum maps obtained at 1mm 
and 3mm with the PdBI show strong emission in 4C 31.04 
(Fig.[T^-b). The continuum distribution at 1mm is fitted by a 
~40 mJy point source centered on the AGN and a Gaussian 
source of ~20 mJy flux density and 1 . 1 kpc x 0. 1 kpc decon- 
volved size (suggestive of a disk). The ~160mJy contin- 
uum source detected at 3mm is point-like. The point source 
detected at 1mm and at 3mm is the higher frequency coun- 
terpart of the non-thermal emission detected at cm-wave- 



lengths. The spectral index (a) of the total continuum emis- 
sion between 84 GHz and 218 GHz (with Sj^^u~^) con- 
firms that the emission is dominated by non-thermal pro- 
cesses: a>l (see Fig.[Tl)). However, the extended disk de- 
tected at 1mm betrays thermal emission from dust. In partic- 
ular, the 1mm disk is tightly linked to the disk identified in 
the R-H color image of Perlman et al. (2001) (Fig.[Tt). The 
dust mass budget (Mdust) of the disk of 4C 31.04 has been 
derived by GB07 by fitting the fluxes measured by IRAS 
(at 100/im and 60/im) and the PdBI (at 1mm). The best fit 
gives a value for M^^^^^ of (4-7)x lO^M©. With a standard 
molecular gas/dust ratio ~100, GB07 estimate that the total 
molecular gas mass of the disk is M^as~(4-7)x lO^^M©. 
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The existence of a massive gas reservoir in 4C 31.04 is 
confirmed by the detection of significantly strong HC0+(1- 
0) emission coming from a rotating molecular disk (Fig.[T]:). 
The HCO+(1-0) disk, of--1.4kpc-size, is detected both in 
emission and absorption. As shown in Fig.[T]: and Fig.[2l 
HCO+(1-0) is detected in emission over ~950kms~^ in 
two separate velocity intervals: channels I and III (as de- 
fined in Fig.O. The strongest emission (III) peaks NE of 
the AGN (offset A in Fig.[T]:). Emission in I peaks SW of 
the AGN (offset B in Fig.[T]:). In contrast, the line is de- 
tected in absorption over ~250kms~^ against the contin- 
uum source of the AGN in channel II (as defined in Fig.O. 
Assuming that the total width of the HCO+ emission pro- 
file is twice the projected rotational velocity (Vrot), GB07 
derive Vrot ~500kms~^ for an edge-on disk geometry and 
a value for the systemic velocity (Ysys) that corresponds 
to z=0.0602. The gas seen in absorption is globally blue- 
shifted by ~ 150 km s~^ with respect to the value of Vsys de- 
rived from the HCO+ emission profile. This result indicates 
that the gas causing the absorption is subject to strong non- 
circular motions. Interestingly, GB07 report on further ev- 
idences of morphological distortions identified in the disk. 
The SW side of the molecular disk is offset from the dust 
disk, a result suggestive of a tilted or warped distribution 
for the gas. 

The HCO+ line emission is a good tracer of the dense 
molecular gas content (n(H2)> lO^cm"^). GB07 derive that 
the mass of the dense molecular gas in 4C 31.04 amounts 
to ^4xlO^M0. In particular, GB07 estimate that the to- 
tal molecular gas mass would range between ~5x lO^M© 
and ~4xlO^^M0, a result consistent with the value de- 
rived from the continuum emission of dust. Most remark- 
ably, the molecular gas masses derived from mm line and 
continuum observations of 4C 31.04 are one-to-two orders 
of magnitude higher than those inferred from HI absorption 
line studies or from optical extinction measurements in this 
CSO (Conway 1996; Perlman et al 2001). As such, this fla- 
grant disagreement highlights the need of using mm inter- 
ferometers to unbiasedly probe the total gas content of radio 
galaxies. 



3 Evolved radio galaxies: the case of 3C 293 

3C 293 is a nearby (z~0.045) moderately large (~200 kpc) 
FR II radio galaxy located at the nucleus of the elliptical 
vv5-33-12. The inner ~4 kpc structure of its radio source 
has been resolved into a jet feature that runs mostly east- 
west (e.g., Beswick et al. 2004). Extensive dust lanes are 
identified in the HST images (Martel et al. 1999). The mod- 
erate (~3'0 resolution CO interferometer map of Evans et 
al. (1999) shows that molecular gas is distributed in a hardly 
resolved disk around the AGN. Using the WSRT array, Mor- 
ganti et al. (2003) report the detection of a broad (^1400 
kms~^) HI absorption feature. The absorption profile iden- 
tifies a surprisingly fast outflow of neutral gas in 3C 293. 




Fig. 3 The 3mm-continuum PdBI map of 3C 293 (con- 
tours: 20 to 200 in steps of 20 mJy), showing the AGN 
source and the jet component to the east. See GB08 for de- 
tails. 



We have recently used the PdBI to probe the gas con- 
tent of 3C 293, observing the HCO+(1-0), HCN+(l-0) and 
^^C0(2-l) lines and their continuum emission at spa- 
tial resolution (Garcia-Burillo et al. 2008, hereafter GB08). 
We have detected the continuum emission at 3mm and 1mm 
not only from the AGN core, but also from the hot spot 
component of the inner jet, located ~1.2'' east of the core 
(Figs. [3] and 131). This hot spot coincides with the contin- 
uum source 'El ' detected at cm wavelenghts by Beswick 
et al. (2004). The molecular gas disk feeding the AGN has 
been mapped and fully resolved in the C0(2-l) line, seen 
in emission everywhere in the disk but towards the AGN 
source (Fig. The disk, of ~7 kpc diameter and signif- 
icantly high mass content (~5x lO^M©), is roughly sym- 
metric and displays a very regular rotating pattern around 
the AGN (see GB08 for a detailed description). 

We have also detected a complex system of HCO+ and 
HCN lines seen in absorption against the AGN mm-core 
and against the eastern hot spot component of the radio 
jet (Fig. [5]). There is a narrow absorption line towards the 
AGN (absorber II), superposed on a much deeper and sig- 
nificantly broader absorption profile (absorber I), as shown 
in Fig. [21 We have also detected a very narrow (<20kms~^) 
absorption system towards the hot spot component (called 
absorber III in Fig. [Sj. These narrow absorption lines are 
likely produced by molecular clouds which are located at 
large radii in the disk. In this scenario, the clouds are pass- 
ing in front of the jet, but they are seemingly not interacting 
with the radio plasma. In this context, it is worth noting that, 
down to the limit of our sensitivity, we have not detected the 
molecular counterpart of the fast HI outflow. 

4 Discussion 

The mm line and continuum observations of 4C 31.04 and 
3C 293 have revealed the existence of significant amounts 
of molecular gas in the two radio galaxies. The derived gas 
masses are comparable to those identified in ULIRG (San- 
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Fig. 4 We overlay the Imm-continuum PdBI map of 
3C293 (contour levels: 10, 15, 20 to 90 in steps of 10% 
of the peak absorption —110 mJy) with the ^^00(2-1) 
integrated intensity emission map (grey scale shown in 
Jy km s~^ units). The eastern component of the jet is also 
detected at 1mm. See GB08 for details. 



ders et al. 1991). Although of similar gas mass, the molecu- 
lar disks of 4C 31.04 and 3C 293 show significantly differ- 
ent properties. While drawing firm conclusions is still pre- 
mature, it is tempting to link these differences with the fact 
that these two sources are thought to represent two different 
evolutionary stages of the radio-loud phenomenon. 

The comparatively smaller size of the molecular disk 
of 4C 31.04, implies that the surface density is an order 
of magnitude higher in this CSO compared to 3C 293. In 
this context, it is worth noting that the C0(2-l) emission 
is barely detected in the PdBI map of 4C 31.04. This is in 
clear constrast with the result obtained by GB08 in 3C 293. 
The previous PdBI observations at 1mm published by GB07 
suggest that the HCO+(1-0)/CO(2-1) is surprisingly large 
in 4C 31.04. Assuming a CO(2-1)/CO(1-0) brightness tem- 
perature ratio of 1, we infer that the HCO+(1-0)/CO(1-0) 
ratio in 4C 31.04 would be ~0.3, nearly an order of magni- 
tude larger than the ratio measured in 3C 293. This indicates 
that the fraction of dense molecular gas is extremely high 
in 4C 31.04 compared to a normal galaxy, and seemingly 
larger than that estimated for 3C 293. This is consistent with 
the comparatively higher gas surface density measured in 
4C 31.04, as derived from the molecular disk sizes of these 
radio galaxies. 

The analytical models of Carvalho (1998) predict that 
a clumpy medium with M^as^^lO^-lO^^M© is able to con- 
fine a radio source on the scale of 0.5-1 kpc. In the case 
of 4C 31.04, the bulk of Mgas is seen to be in a disk per- 
pendicular to the radio source's axis, a geometry that would 
decrease the efficiency of frustration, however. Within the 
frustration scheme, it is foreseen that a disturbed ISM would 
reflect the interaction between the radio plasma and the den- 
se confining cocoon. In the case of 4C 31.04, there is evi- 
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Fig. 5 We show in this figure the HCO+(1-0) p-v plot 
taken along the apparent axis of the jet at 3mm (roughly 
oriented east- west). We highlight the presence of a narrow 
(<30kms~^) absorption system on the east side of the jet 
(absorber III). There is also a narrow absorption line at 
v=250kms~^ that betrays absorption against the AGN, su- 
perposed on a much deeper and broader absorption detected 
around v=60kms~^. Levels are 4, 6, 8, 10, 20 to 80 in steps 
of 15% of the peak absorption —180 mJy. See GB08 for 
details. 



dence that a radio lobe-ISM interaction may be at work. The 
HST image of Perlman et al (2001) shows cone-like fea- 
tures aligned with the radio axis of the source (see Fig.[T]:), 
the likely signature of gas shocked by the jet. Moreover, 
the distribution and kinematics of molecular gas probed by 
the HCO+ observations discussed by GB07 illustrate that 
the detected rotating disk is not in a fully relaxed state. The 
revealed distortions indicate that the disk is still settling af- 
ter a merger or an event of gas accretion. Alternatively, the 
jet and the cone-like features may be interacting with the 
disk, thus producing the reported distortions. In the case of 
3C 293, the distribution and kinematics of molecular gas are 
reminiscent of a regularly rotating and dynamically relaxed 
disk. Contrary to the case of 4C 31.04, there are few signa- 
tures, if any, of an ongoing interaction between the jet and 
the molecular gas disk in 3C 293. 
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